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Abstract. In the [BCS] paper the theory of superconductivity was developed for the BCS model, 
in which the (instantaneous) interaction is only between fermions of opposite momentum and spin. 
Such model was analyzed by variational methods, finding that a superconducting behavior is ener- 
getically favorable. Subsequently it was claimed that in the thermodynamic limit the BCS model is 
equivalent to the (exactly solvable) quadratic mean field BCS model; a rigorous proof of this claim 
is however still lacking. In this paper we consider the BCS model with a finite range temporal 
interaction, and we prove rigorously its equivalence with the mean field BCS model in the thermod- 
inamic limit if the range is long enough, by a (uniformly convergent) perturbation expansion about 
mean field theory. 



1. Introduction and main results 

Bardeen, Cooper and Schreifer [BCS] developed their theory describing superconductors by the 
BCS model, in which the interaction has infinite range and only fermions of opposite momentum 
and spin (Cooper pairs) interact; the Hamiltonian of this model is 

H B cs = J2j v d ™t (- H) a t - j v **>t, a i- J£ J v d yvi-A 

where CT are creation or annihilation fermionic field operators with spin a in a d-dimensional box 
with side L and V = L d , m is the mass and A > is the (attractive) coupling. By a variational 
procedure it was found that it was energetically favorable to form a superconducting phase. Later 
on it was realized that the properties of such superconducting phase are identical to the ones of 
the mean field BCS model, an exactly solvable model in which the interaction is quadratic and the 
Hamiltonian has the form 

H MF = I A| 2 + E / y d ™h (-||) 4. - VAA[E J y dya^a^J ^A[^ jf dya± a a~_ a ] 

(1-2) 

where A is a complex number to be determined minimizing the ground state energy (that is A 
solves the BCS gap equation). It has been argued in several papers, starting from [BR],[B],[H],that 
in the limit V — > oo the reduced BCS model (1.1) and the mean field model (1.2) have the same 
correlation functions; this seems quite natural also by analogy with lattice classical statistical 
mechanics in which infinite range interaction gives mean field behavior in the thermodynamic 
limit. Indeed many arguments has been given to support this claim in the last fifty years but, as 
far as I known, a rigorous proof is still lacking; aim of this paper is show that a simple proof of 
this claim can be given at least if the instantaneous interaction in the reduced BCS model (1.1) is 
replaced with a long (but finite) range time interaction 

We consider then a generalization of the reduced BCS model in which fermions are on on a cubic 
lattice with step 1 and a time- dependent interaction between Cooper pairs is considered; indeed, as 



stressed for instance in [CEKO], a realistic model for superconductivity should include a bosonic 
Hamiltonian describing phonons and a boson-fermion interaction, which can be written in a purely 
fermionic model only if time dependent interaction between fermions is included. The two point 
Schwinger function of the reduced BCS model on a lattice with time dependent interaction can be 
written as Grassmann functional integral in the following way 

/P(d^)e- V - h ^.«J d ** s -^ ; -- 
where J dx = J 2 ^ dxo J2xeA ano - A is a d-dimensional lattice with step 1 and 
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(1.4) 



In the above expression {V'ker} i s a se t of Grassmannian variables, k € T)l$ where £>l,/3 = T^l x ^ p, 
with V L = {k = 2im/L,n e Z d ,-[L/2] < n t < [(L - l)/2]} and Vp = {k = 2(n + 1/2)tt//3, n e 
Z, — M < n < M — 1}, and P(dip) is a linear functional on the generated Grassmann algebra such 
that 

/ n#)C llCTl V£ 2 , CT2 = £ d /3<5k llk2 <W 2 2(ki) , 5(k) = _ lkQ+ \ {k) _^ ■ (1-5) 

where 

d 

e(fc) = ^(1 -cosfci) (1.6) 

is the dispersion relation and /i is the chemical potential. We define also Grassmannian Held ip^ 
is defined by, if x = (a;, x ) with x G (— f , f ] and a; = (x\, ..,x d ) with £j = 1, 2, .., L, 

<- = ^9 £ ^ e±4k ' X (1-7) 

The external field ft, is introduced to break the number symmetry and which will be removed after 
the thermodynamic limit L — > oo will be taken and v(x — y ) is a Kac potential with a long but 
finite range potential for definiteness we choose 

n =0,±l,±2,..,±M 

Finally M is an ultraviolet cutoff in the time direction introduced to make the Grassmann integral 
well defined, and and the limit M — > oo must be taken before the thermodynamic limit V — > oo. 
As we mentioned above, it was claimed in [BR] that the reduced BCS model (not solvable) should 
be equivalent (in the sense that the Schwinger functions coincide) to the mean field BCS model 
(solvable) in the limit L — ► oo. In [BZT] indeed it was shown at each order of the perturbative 
expansion that the difference of the correlation functions between the reduced and the BCS model 
goes as 0{V~ l ) at each order, but to make this argument rigorous one has to prove the uniformity 
of the convergence of the perturbative expansion. A similar perturbative argument in a more 
modern (RG) language has been given in [SHML], in which it is pointed out the similarity of 
the perturbative expansion of the reduced BCS model with the so called O(jj) expansion. In 
[B] and [H] the proof of such equivalence was based on the idea that the spatial averages of field 
operators like V~ l J dxa^ s a^ a - may be substituted by numbers in the thermodynamic limit, since 
commutators with them has an extra one volume factor. Such a proof was criticized by several 
authors; for instance in [TW] it was shown that the convergence of the reduced BCS to the mean 
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field model is true only in a rather small subspace and not in general. In [M] a new proof of 
the equivalence based on a functional integral approach was given, but in the analysis involves 
unjustified exchange of the L — > oo limit with the M — ► oo limit. Finally in [T] a correct proof of 
such equivalence was given, but only under that rather unrealistic assumption the the dispersion 
relation is a constant (degenerate BCS model). 

It is apparently surprising the difficulty in proving that a infinite range interaction interaction like 
the one in (1.1) leads to a mean field behavior in the thermodinamic limit; indeed in classical statis- 
tical mechanics for spin lattice systems the proof of a similar statement is a two line computation. 
The difficulty in the quantum case can be clearly understood in the functional integral formulation 
(1.3); in such a representation the interaction V is not factorized contrary to what happen in the 
Hamiltonian formulation, and this make the model not exactly solvable. Of course by replacing 
v(xo — yo) in (1-3) with a constant (that is we consider an infinite long range time interaction 
= oo) the interaction in the functional integral is factorized and the model is exactly solvable; 
mean field behavior in the thermodinamic limit is then easily established, by performing a saddle 
point analysis essentially identical to the one for long range spin systems, see [L]. 

Aim of this paper is to prove that even if the range in (1.8) is finite, so that the interaction 
is not factorized and the model not solvable, the BCS model (1.3) is equivalent to the mean field 
BCS model if n is large enough, in the limit V — ► oo; that is the BCS model has a phase transition 
into a superconducting state described by the BCS theory. 

Our main result is the following. 

Theorem Assume fi < 2 and A > 0; there exist (3 C (\) and n ((3) > such that for (3 > (3 C (X) and 
< k < kq([3) the Schwinger functions (1.3) with interaction (1.8) are such that 

i — ik + e(k) — a , 
lim lim < ^Zj>t„ >= ~ — (1-9) 

k,CT fcg + ( £ (fc) - M ) 2 + A|A| 2 

lim lim < VwV^w >= _VAA ^ ^ 

fc-»0+L^oo ' ' k 2 + (e(jfe) - /i) 2 + A|A| 2 

where A = A(/3) is the positive solution of the BCS gap equation 

^ r d k tangh(I^E 2 (k) + XA 2 ) 
J i^Y 2((e(fc) - fi) 2 + AA 2 ) 

and (3 C (\) is the minimal (3 for which (1.11) admits a solution. 



The above Theorem ensures that, at a fixed temperature (3 and for range large enough, the 
BCS model has the same behavior of the BCS mean field model; in particular for [3 > [3 C a gap is 
generated and the particle number symmetry is broken as < tp£ a^-k a 1S different from zero; this 
means that there is a phase transition into a superconducting phase for temperatures low enough. 
As an immediate corollary, it follows that for an interaction like (1.8) with an exponentially large 
range O(e^) a gap is generated and the particle number symmetry is broken for temperatures 
small enough. 

The proof of the above statement is by perturbation theory about the mean field theory, using 
as a perturbative parameter the inverse range k of the Kac potential (1.8); this is a a classical 
approach in classical statistical mechanics to prove phase transition beyond mean field theory, see 
for instance [LMP]. We will show that the correction with respect to the mean field Schwinger 
function is expressed by a convergent series expansion (uniformly in the volume) and each order is 
0(V~ l ); uniform convergence is established via determinant bounds for fermions. We can prove 

convergence only for small n, as it turns out that k — 0{[3~^~); of course it would be very 
interesting to prove convergence up to «o = 0(1) or even for any k, so otaining a real solution of 
the BCS model with instantaneous interaction. 
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2. Partial Hubbard-Stratonovich transformation 

In momentum space we can write the interaction V in the following way 

V = ~T^¥^ ^ $>(Po)< k ^,_ k+Po Vv,k'C CT ' ; _ k ' +Po (2.1) 

VH > k,k' P0= 2*L ( „ + i) a,a> 

where we have used that p = fc ,i + fco,2 = ^f-((«o,i + ^0,2) + !)• We split the interaction V as 
sum over two terms 

V = V + V 2 (2.2) 



^ = ~WfY S I>W<k^,-k+ Po ^,k^- CT ',-k' +Po (2-4) 

^ ^ k,k',| P0 |>^CT,<7' 

Note that V can be written as, e = ± 



V = -A+A" A e = 



p£ = ppEE^-M (2-5) 



that is can be written as the product of the total number of Cooper pairs. Let us consider the 
generating function of the Schwinger functions 



e SL,f> 



= J P(d^)e 2A+A "-*e- fc ^ A+ - h ^ A "e/ & S. [J - + .^-- + <' J -- ] (2.6) 



where J are external Grasssmann field, so that 

<^n,>=Q^ F S{J)\. J= , (2.7) 
x y 

By using the identity (Hubbard- Stratanovich transformation) (4> = u + iv, 4> = u — iv, u,v € R) 

e 2ab = — [ dudve- 1 ^ 2 e a * +b * (2.8) 

27T J R 2 V ; 

we can rewrite the above expression as 

e s L ,e, h {J) = 1_ f dudve -m 2 /p(^ )e -v e (0-^)A + +(^^)A-J d x E jj+^-„ + <^-j 
2?r J R 2 J 

(2.9) 

Performing the change of variables (u, v) — > y '/3V (u, v) we obtain 

e SL,0, h {J) = PL f dudve'^ (v2+(u+ ^ )2 e- 0V:FL -f- hiu ' v)+BL -f- hiu ' v '' l>) (2 10) 

27T Jr2 ^ ' ' 



where 



-/3VJ^ L , , h (u,v)+B L ^, h (u,v,J) _ J p(rf 7 /,) e -V e yA0X)+ + yA^X) g J dx [J+ CT ^ X <T +^ + <T J x (2.11) 



and by definition BL,p,h(u,v, J) is vanishing for J = so that TL,p,h{u,v) is given by 

e -0vr L ,p, h (u,v) = J p ^y-v e V\4>v++V\4>v- (2 .12) 

We are interested in computing the two point Scwhinger function, given by (2.7) 

< ^-ee^a >= [ dudve~^ ^ u +^) e -0V^, , hM ge j (k> ^ y) (2 . 13) 

where S L ,p,h{u,v) = djedje'B(J,u,v)\j= and 

Z L0h = [ dudve-^ v2+{u+ Vz )2 e-P v:F ^.^ (2.14) 

We will show in the following section that 

F\,0,h(u, v) = t B cs + FL,p,h{u, v) (2.15) 

where, if E{k) = e(k) — fj, 

is the free energy in the mean field BCS model [BCS] and f~L,p,h is the perturbation to the mean 
field; we will show in the following section that, for < k < Ko(/3), Ko([3) = C~ x (5^^~ , for a 
suitable constant C 

\^L t0< h{u,v)\<C^{K 2 O^ d+2 (2.17) 

hence V^L,0,h(u, v) it is uniformly bounded as V — > oo; it is more convenient to call Vf~L,p,h(u, v) = 
^L,0,h(u,v) and we can write the two point Schwingcr functions as 

— L_ f W !)e -^l" 2 +(''+^) 2 + tscs ( v »l e -^.^( v '5[' £ ( ;(k,«,«) (2.18) 

Z L,/3,h JR 2 

By the saddle point Theorem, for (3 large enough 

-f3L(v 2 + (u+-±=) 2 +t BC s(u,v)) 

lim j— = 6(u)6(v - v ) (2.19) 

f dudve^ L{v+{u+ ^ )+tBCs{u ' v)) 

where Vq is given by the negative (for h > 0) solution of 



(2ir) d 2(E 2 (k) + \v$) 
In the limit h — > it reduces to the BCS equation (1.11). Moreover we will show in the following 



section that S e £ £ g — S^ £ 0' BCS is 0(V 1 ) so that the Theorem follows. 



3. Convergence of series expansion 

3.1 The partition function 
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We can "absorb" the quadratic fermion term in the the free interaction 

J P{di,)e^ v++ ^ v 'e-^ = e - 0vtBCS J P a {d^)e~^ (3.1) 



where 



and 



^WO = ~ E / d * d yv(*° W))^^^, (3-2) 



^o-yo) = lj2e^^- ¥ (3.3) 



(3.4) 



and, if (J — VA0 

w) = n^(-^E E 

where 7V(k) is the normalization of P CT (df/0 and 

i v fc 2 + ff 2 (fc) + H 2 
-^L lo s fc 2 + ^( fc ) (3 - 5) 

and the 2x2 matrix T(k') is given by 



We can equivalently write, see for instance [L], tscs as (2-16) and of course t\ < y|A|C[l + |^|]. 
The propagator of P a (dip) is given by 

k 

We decompose the free propagator g k into a sum of two propagators supported in the regions of 
k "large" and "small", respectively. The regions of ko large and small are defined in terms of a 
smooth support function H (t), t£R, such that 

with 7 > 1. We define h(ko) — H (\ko\) so that we can rewrite <?k as: 

3 k = fl^-)(k)+^(k) (3.9) 

where 

^:;°(x,y) = -^E e " k(x_y) ^ fc «)[ T_1 ( k )]^' ( 3 - 10 ) 



k 

(u.v.) j 



9^ ) (x,y) = ^E e " k(X_y) ( 1 -^o))[T- 1 (k)] £ , e , (3.11) 
In the Appendix we show that 

P a (dr- v -)e-^ ir +^ = e - v °^ ir ) (3.12) 



with 

V ° = -yYl / dxd y i ( x O-yoHt°' l l J t-^y^y,-<J>+Yl / dXl " / ^"^n ( Xl ' ' ' ' X2 «) II ^ ,ct 
cr.<r' n=l i=l 

(3.13) 

with 

^ Jdx 1 ---dx 2n \W^\x 1 ,..,x 2n )\ ^^lAr^ 1 '"- 1 )^ 3 ) 1 "^ 1 '"- 1 ) (3.14) 



3.2 Convergence of the infrared integration 

We define a distance d(x,y) i/3 = (^(^o, J/o), d L (xi, yi), ...,d L (x n ,y n ) as 

d/3(x ,yo) = -sin^(x - y ) d L (xi, yi) = — sin yfa - y{) (3.15) 

7T p TT L 

In order to perform the infrared integration we need the large distances behaviour of the infrared 
propagator. 

Lemma For any integer N the following bounds hold 



\^y)\<P 1 + [ ^%_ y)]N 



InMfxvKfl 9lL (317) 



Proof. The above bounds follows by integrating by parts. Consider integers N , Ni, .., Nd note 
that, i = l, .., d 



dh {x l ,y i ) N% dp (x , yo) °9i!e' (x - y) = 



^i^L^N.+^N^^No+N^^-M^Q^QN, [(1 _ /i(fc ))[T - 1 (k , )]e,e'] , (3 ' 18) 

where dk and dk denote the discrete derivatives. The bound then easily follows noting that 
T ~ 1 (k / )] e e ' is bounded by C(3 and each derivatives over it is bounded by an extra (3. The non 
diagonal term has an extra k 2^ g |2 m the bound, from which we sec that the bound is uniform in 
a.)M 

We can write 

J P(d^)e' v °^ = eEriof- 1 )"^^-^) (3.19) 
where £ T are the fermionic truncated expectations 

l0 § / P(d^)e XXW \ x=q (3.20) 

We write (3.13) as 



[ dx P W(x P )4>(P) (3.21) 

p 

where P v is the set of field labels appearing in (3.13), W(xp^) are the kernels in (3.13), that is 
V^ 1 v(x - j/o ) or W(xi, ..,x 2 „) and 

m = n €\ f fUf) ( 3 - 22 ) 



Then we get 

£ T (V°;n)= J2 [ dx Pi- [ d XPn W(x Pl )...W(x P J£ T (iP(P 1 )....t>(P n )) (3.23) 

The fermionic truncated expectations can be bounded by the formula (see [GM] for example), if 
s > 1, 

£ T ^(P 1 ) 1 ...MP s ))=Y J T{9e l ,s [ {*i~y l ) f dP T (t)dctG T (t) , (3.24) 

where 

m = n </!,.(/) ( 3 - 25 ) 

feP 

and 

a) T is a set of lines forming an anchored tree between the cluster of poins Pi,..,P s i.e. T is a set 
of lines which becomes a tree if one identifies all the points in the same clusters. 

c) t = {tiy e [0, 1], 1 < i, i' < s}, dP T (t) is a probability measure with support on a set of t such 
that ti.i> = Ui ■ Ui> for some family of vectors e M s of unit norm. 

d) G T (t) is a (N - s + 1) x (N - s + 1) matrix, 2N = \Pi\ + .. + \P S \ whose elements are given 
by G fj,i>j> = ti,i>9 u -, u +(xij - Yi'j') with {fijj^,) not belonging to T. 

If s — 1 the sum over T is empty, but we can still use the above equation by interpreting the r.h.s. 
as 1 if Pi is empty, and detG(Pi) otherwise. 

We bound the determinant using the well known Gram-Hadamard inequality, stating that, if M 
is a square matrix with elements My of the form My =< Ai,Bj >, where Ai, Bj are vectors in a 
Hilbert space with scalar product < ■, ■ >, then 

|detM| < JIPill-H^H . (3.26) 

i 

where || • || is the norm induced by the scalar product. 

Let H = M. s ® Hq, where Ho is the Hilbert space of complex four dimensional vectors F(k) = 
(i<i(k), . . . , F 4 (k)), Fi(k) being a function on the set 2?-,-, with scalar product 

< F, G >= J2 i- R E F * ( k ) G *( k ) • ( 3 - 27 ) 



r ^ k 



and one checks that 

?t _+. ..„(x) 



G ij,i'j> = - Yi'j') =< ^ ® ^4 x(/ - );W(/ - ) ,u i , ® B x(/ + ,) > , (3-28) 

where e R s , i = 1, . . . , s, are the vectors such that ti^ — • u^, and 



A x u (k) = c tk ' x VMM . / (-<*o + 0, fj, 0), if w = +1, 



y^+^TFF 1(0^,0,1), if w = -i, 

R = p ik' y a/w r (i,i,o,o), ifw = +i, 

v^q-^TM^' 1(0,0,1,-^0-^))), if w = -i. 

Hence from (3.26), as < C and < we find 

|GS W l<^- s+1 ^-' +1 (3-30) 



where C\ is an 0(1) constant. 
By using the above formula in (3.23) we get 

|£ T (V°,..,V°)|< ]T f dxPl ... /dxpJ^(x Pl )|...|^(xpJ|^[J]irW(x ; -y i )l] 

Pl,...,P„ J J rp ^j, 

(3.31) 

where we have used that / dPp(t) = 1. The number of addends in ^ T is bounded by nICJ ■ 
In order to bound the integration over propagators we use antipcriodicity 



1-0/2 r- r0 

dx g s ,e'(x,Xo) = dx g £ ^(x,x ) + dx g £ ^ (x, x ) = / dx Q (g £t£ ,(x 1 x a )+g £ , £ ,(x ~(3) 

J-0/2 J\x o \>0/2 J-0 



The tree T realizes a connection between all the V, and we get the bound 
r n 1 

/ n^Etni^w^-y)!] ^ w)p (n - i){d+i) (3.32) 

In order to perform the integration over the remaining coordinates we note that if W(xp) = 
V~ 1 v(x ) then 

\Hxo)\ = A E ^^^i^E^^-w ( 3 - 33 ) 



' <=o^° u ^ n#0 



so that 



y dxV _1 |i;(a;o)| < C(k(3) 2 (3.34) 

if Cis a suitable constant. On the other hand if W(xp) — W(xi, ..,x 2 „) we use the bound (3.14); 
then we get, assuming k/3 < C~ x in order to sum over Pi 

n 

i=l Pi 

{l3V)n\C n X n {K 2 l3 3 ) n f3 {d+3)n l3- {d+1) < (f3V)C\(K 2 (3 d+6 ) n /H d+1 )n! 
Finally the following bound can be found, calling Th,p,h — tBCS + Fhfi,h 

\T L , fj .h\ < CA( K 2 /? d+6 )/T (d+1) 

assuming that k < C" 1 (3^~ = kq(/3) to assuring the convergence of the sum over n. 

Remark The above analysis immediately imply a bound for the effective potential J P(dip^^^ )e~ v ^ + ^ 
where <j> is an external fermionic field. The kernels of the effective potential W^ n )(xi, ..,x„ e ) at 
order n obey to the bound 

^ J dxi..dx r ,|w n (x 1 ,.., Xn )i < c n \ n ( K 2 (3 d+6 rr^r id+1) (3.35) 

as now the propagators are 2n — ^ 
3.3 Extracting a volume factor 

The above analysis says that FL,p,h, which is the correction to the mean field, is given by a 
convergent expansion for sufficiently long range interaction < k < Ko(/3). We prove now that we 
can improve the above bound by a factor V^ 1 . 
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Consider first the case in which in (3.23) there is at least a W(Pi) associated to v. We can write, 
by using that for the fields in £ T holds the rule ip x = j dx! g(x — x') 

-^-^£ T (V°;n) = —.^tj{y7 [ dx 1 dy 1 dx a dx b dx c dx d v(x ^ 1 - y .i) 
n\ pV n\ pV V J 

9+,e a (xi -x a ).g +i£b (xi -x 6 )g_ i£c (yi - x c ) ff _ :£d (yi - x d )H^ b ^ £d {x a , x 6 , x c , x d ) 
J dyLidy 1 dyL a dyL b g + -(xi-yi)v(x aA -yoA)g + ,c a (xi-Xa)g+,eAyi- x b)H^ (3.36) 

^ V j d * ld y idXadXb 9-c ,-e(0)v(X0,l - 2/04)9£',£a( x l - x a).9e',£(,(yi - x b)H 2 ^ n £b (x a , X b )} 



e=± 

where 



^L,e d (x a ,x 6 ,x C) x d ) = ^^^^ ( y;n-l) (3.37) 

«(x ,x b) = ^^ ( F;n-l) (3.38) 

Note that the last addend in (3.36) (corresponding to a tadpole contribution) is vanishing; in fact 
it can be written in momentum space as 

E V,o^o)^£^(k')^, £6 (k' +p )H^ Eb (k>, Po ) = 

po/0 k' 

The first addend in (3.36) can be bounded in the following way, remembering that |<7 £>£ /(x, y)| < (3 
C^/? _1 (k;/?) 2 /? 2 sup[^ dxi| ff+!£a (xi -x a )\]sup[J rfyi.g-, £c (yi - x c )] 



< 



1 f 

n\~]iV 



j dx a dx b dx c dx d \H^"l ec Ed (x a , x 6 , x c , x d )| (3.39) 



The bound for the last integral is given by (3.35) with n e = 4; hence the first addend in (3.36) 
obeys to the following bound 

C"^ r /3- 1 ( K /3) 2 /3 4 /3 2(d+1) (« 2 /3 d+6 )"" 1 r 2 r (d+1) (3-40) 

so that summing over n wc have the bound v (np) 2 f3 d + 2 . 
Finally the second addend in (3.36) can be bounded by 



(k/3) 2 /3sup[ / dxi|5 +)£o (xi -x )|]sup[ / dyig+, eb (yi - *b)] / dx a dx b \H 2 ^ n E b (x a ,x b )| (3.41) 

x a J ii J J 



and again using (3.35) with n e = 2 we get that hence the second addend in (3.36) obeys to the 
following bound 

A™ 

C n —f3- 1 (^) 2 ^f3 2id+1) (K 2 p d+6 ) n - 1 r 1 r {d+1) (3.42) 

Of course there is no v, we can apply the same reasoning to one of the kernel W(xi, ..x„). Summing 
over n we have the bound y (k 2 f3 3 ) (3 d+2 ; then, for k < fc we get the better bound 

\^L,f3, h \<Cy (K 2 p 3 )/3 d+2 (3.43) 
10 



3.4 The integration of S 

By performing the change of variables, if ip = and g is the matrix propagator of P a (dtp), 

V'k — * V'k + gi>k, we get for two point Schwinger function the formula 

S £ , e ,(k,u,v) = ffe>e ,(k) + ^ fle, £ "(k)y 2;£ ,, £ -(k) 5 ^^(k) (3.44) 

s",s"' 

where V 2 - £ i e »< (k) is the kernel of the effective potential with two external fields; it can be bounded 

by 

|*W(k)| <^JdxJ dy\V 2 , £ , £ ,(x,y)\ (3.45) 

By using (3.35) we get that 

|^W( k )l < C n X n (k 2 /3 d+6 ) n /3' 1 /3~ {d+1) (3.46) 
By can improve the above bound as described in §3.3. If there is at least a W = V~ 1; v, we get 

V2 ;e ,e'(x,y) = ^ J dzdz'dx 1 dx 2 dx 3 dx4g(z-x 1 )g(z-x 2 )g(z'-x 3 )g(z-x 4 ,)V 6 (x,y,x 1 ,x 2 ,x 3 ,x 4 ,)+ 

^ / dx J dydx x dx 2 J dzdz' g{z - z')g(z - x 1 )g(z' - x 2 )V^(x,y,xi,x 2 ) 

where we have used that the tadpoles contributions is vanishing. The the integral over x, y times 
(/3F) _1 of the first addend can be bounded, using also that sup\g\ < (3 

J dx J dydxidx 2 dx 3 dx i J dzdz' 

\g(z - x{)g{z - x 2 )g(z' - x 3 )g(z - x 4 )||S(x,y,xi,X2,X3,x 4 )| < 

-^(K(3f(3 2 (sup xi \ J dzg(z-xi)\)(sup X3 \ J dz' g(z' -x 3 )|)-^ J dx J dydx 1 dx 2 dx 3 dx 4 \S{x,y,x 1 ,x 2 ,x 3 ,x 4 )\ 

(3.47) 

By using (3.35) we get for (3.44) the bound y(nP) 2 (3 d+1 . On the other hand the second addend 
is bounded by 

^^P{sup Xl \ J dzg(z - xi)|)(sitp X3 | J dz'g(z' - x 3 )|) J dx J dydxidx 2 |V2(x, y, xi, x 2 )| 

from again by using (3.35) we get the same bound y(n/3) 2 [3 d+1 . If there are no vertex v, we can 
repaet the above argument on the kernels of (3.14). Hence from (3.44) 

|S c , e ,(k,«, V )-<fe, e ,(k)| < ^( K 2 /? 3 )/3 d+1 (3.48) 

Appendix Al. The ultraviolet integration 

The integration of the ultraviolet part (3.12) can be done by a multiscale analyses; it is quite 
standard and we refer to §3 of [BM] for details in a similar case. It is convenient to introduce an 
ultraviolet cut-off N by writing 

N 

^](x) = ^ 9 ( fe )(x) (3.1) 
k=i 

where 
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with h k {ko) = H (j- k \k a \) - Ho{-y- k+1 \k \). Note that limjv-ooS 11 '^*) = g( u - v -\x) and that, 
for any integer K > 0, g^(x) satisfies the bound, for any integer K 



| 5 ^>(x)| < . , CK t (3.3) 

I + (~f k \x \ + \x\) K V ^ 

We associate to any propagator (x) a Grassmann field ip^ and a Gaussian integration P(dil)^) 
with propagator </ fe )(x). We can rewrite as: 

V (o) (0)+F/3^i = - lim log / P(#(°))P(#( 1 ))---P(#( JV ))e-^ [1 ' W1+ ^ (3.4) 

AT— >00 J 

We can integrate iteratively the fields on scale iV, iV — 1, . . . , fc+ 1 and after each integration, using 
iteratively an identity like (3.19), we can rewrite the r.h.s. of (3.4) in terms of a new effective 
potential V [1 ' k] : 

N 

(3.4) =lim o {v/3 E o~ l °S P» (1 ' »)P(# (y) )-P(# (1 ' t) )e- v "''"» [U1 ^} (3.5) 

with V^'^^V'' 1 '^) admitting a representation in terms of trees defined in the following way: 

1) Let us consider the family of all trees which can be constructed by joining a point r, the root, 
with an ordered set of n > 1 points, the endpoints of the unlabeled tree, so that r is not a branching 
point, n will be called the order of the unlabeled tree and the branching points will be called the 
non trivial vertices. The unlabeled trees are partially ordered from the root to the endpoints in 
the natural way; we shall use the symbol < to denote the partial order. 

Two unlabeled trees are identified if they can be superposed by a suitable continuous deformation, 
so that the endpoints with the same index coincide. It is then easy to see that the number of 
unlabeled trees with n end-points is bounded by 4". 

We shall consider also the labeled trees (to be called simply trees in the following) ; they are defined 
by associating some labels with the unlabeled trees, as explained in the following items. 

2) We associate a label k > with the root and we denote T(k,N),n the corresponding set of labeled 
trees with n endpoints. Moreover, we introduce a family of vertical lines, labeled by an an integer 
taking values in [k, N], and we represent any tree t e T<k,N),n so that, if v is an endpoint or a non 
trivial vertex, it is contained in a vertical line with index h v > k, to be called the scale of v, while 
the root is on the line with index k. There is the constraint that, if v is an endpoint, h v = N + 1. 

The tree will intersect in general the vertical lines in set of points different from the root, the 
endpoints and the non trivial vertices; these points will be called trivial vertices. The set of the 
vertices of t will be the union of the endpoints, the trivial vertices and the non trivial vertices. 
Note that, if V\ and v 2 are two vertices and V\ < v 2 , then h vi < h V2 . 

Moreover, there is only one vertex immediately following the root, which will be denoted vo and 
can not be an endpoint; its scale is k + 1. 

3) With each endpoint v of scale h v = N + 1 we associate V (3.2). Given a vertex v, which is not 
an endpoint, x„ will denote the family of all space-time points associated with one of the endpoints 
following v. 

4) The trees containing only the root and an endpoint of scale k + 1 will be called the trivial trees. 

5) We introduce a, field label f to distinguish the field variables appearing in the terms V associated 
with the endpoints. The set of field labels associated with the endpoint v will be called I v . 
Analogously, if v is not an endpoint, we shall call the set of field labels associated with the 
endpoints following the vertex v; x(/), and cr(f) will denote the space-time point, the a index and 
the u index, respectively, of the field variable with label /. 

The effective potential can be written in the following way: 

oo 

VW(^) + L(3E h+1 = J2 E V {h \r,^) (3.6), 

"=1 reT^h,N)n 
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where, if vo is the first vertex of r and t\, ..,t s (s = s Vo ) are the subtrees of r with root vq, 
V^(T,tp^- h ^) is defined inductively by the relation 

and y( /l+1 )(r i , if is trivial and h = N - 1, it is equal to V (2.4). 

By iterating (3.7) we can write V^ h \r, ip(- h>> ) in the following way. We associate with any vertex v 
of the tree a subset P v of I v , the external fields of w. These subsets must satisfy various constraints. 
First of all, if v is not an endpoint and vi, . . . ,v Sv are the vertices immediately following it, then 
P v C UiP Vi ; if v is an endpoint, P v = I v . We shall denote Q Vi the intersection of P v and P Vi ; this 
definition implies that P v = DiQ Vi . The subsets P Vi \Q Vi , whose union will be made, by definition, 
of the internal fields of v, have to be non empty, if s v > 1. 

Given r G T(h,N)n, there are many possible choices of the subsets P v , v G r, compatible with all 
the constraints; we shall denote V T the family of all these choices and P the elements of V T . Then 
we can write 

vW(T,v ( - fc) )= 51 ^ (h) (^P); (3-8) 

y( /l )(r, P) can be represented as 

f< h >(r,P) = | rfx t , V;^' l )(P„ )4 h p +1) (x„ ) , (3.9) 

with K^p^fava) defined inductively (recall that h Vo = h + 1) by the equation, valid for any v G r 
which is not an endpoint, 



<p ) (x,)=^ T n[^ +i )(x, i)] - 

St " i=i (3-10) 



where iv v is an endpoint Ki hv \x v ) is the kernel v(x). We call x _v ertices the vertices v of r such 
that £"£ is not trivial. 

By using the representation of the truncated expectation analogous to (3.20) and the Gram 
inequality we get that the contribution from a tree r G T[i,k\,n associated to a kernel with 21 
external legs can be bounded as (see §3.14 [BM] for details in a similar case): 

-!- f dx 1 ■ ■■dx 2l \W^' N) (T;x 1 ,a 1 ,e 1 ; . . . ; x 2i , a 2U e 2 i)\ < 

" "" (3.H) 



< C"|A(K/3) 2 r7- fe ( n - 1 ) Yl 1 



{h v -h v ,)(n v -l) 
v not c.p. 



where v' is the vertex immediately preceding v on t, n v is the number of endpoints following v on 

T. 

Note that the "dimension" n v — 1 is vanishing when when n v = 1, so that the above bound is not 
suitable to sum over trees. However in the case n v = 1 the above bound can be improved. If fact 
let be v a vertex with n v = 1 such that its set of internal lines is not empty; to such vertex is then 
associated a truncated expectation 

£l (V) = ~ J dx J ^(xo - yo)O y , CT ^+-(x, y) (3.12) 



a, a 1 
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where we have used that the contration of or t()~ij)~ in V is vanishing by momentum 

conservation (remember that po = 0). In (3.11) <?+l(x, y) is bounded with a constant; however 
such bound can be improved by writing 

9+L (x;y) = g+L(x,x ; y,x ) + (x - y ) / dtd xg g^L(x,x ;y,yo + t(x - y ))) (3.13) 

Jo 

and noting that 

l*-(x,x)| < |Jg Ze-'l'AM ki + E ^ ) + - \ + E ^ - || + H , I < C7 

(3.14) 

Then the contribution to the kernel of (3.12) coming from the first addend in the r.h.s. of (3.3) is 
bounded by C7 - ' 1 " ; in the contribution coming from the second addend we bound the propagator 
by C7"' 11 ' so getting a factor V -1 J dr|ro||u(ro)| < C(3, so that at the hand we get a bound for 
(3.12) C7? 7 -V 
Then we get the bound 



-h v 



— J dxi ■ ■ ■ dx 2 i\W^ ,N] (T;x 1 ,(T 1 ,e 1 ; . . . ; x 2i , <t 2 i, e 2 i)\ < 



<C"|A( K /3) 2 r 7 - fc ("- 1+ " tad ^" \{ 7 -(^-V)(",-i+*„) 5 

v not e.p. 

where z v = 1 if n v = 1 and otherwise, and n tad is the total number of c-vertices v with n v = 1 
such that its set of internal lines is not empty. Then, proceeding as in §3.14 of [BM], one can sum 
over t and the bound (3.14) is proved. 

Finally by proceedig as in §3.3 it is easy to see that we can extract a factor 0(y _1 ) from each 
kernel W. 
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